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Abstract. We show that the structure of the quantum group su q {2) is in- 
timately related to the random walks on a two dimentional lattice. Using this 
connection we obtain an appropriate candidate for the exact area distribution 
of closed random walks of length N on a two dimensional square lattice. We 
compare our results with exact enumeration. 
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1. Introduction. Let us consider a spinless electron on a two dimensional 
lattice and submitted to a uniform magnetic field along the z-direction and 
perpendicular to the plane of motion. The Hamiltonian of the system is: 

h=±(? + eAf. (1) 

The system is not invariant under translations but there is an invariance under 
the so-called magnetic translation operators [1]: 

w(a) = exp(-d.k), (2) 

where a is an arbitrary two dimensional lattice vector and k is: 

k = k x x + k y y = p + eA + er x B. (3) 

There is a relationship between magnetic translation operators, quantum groups, 
Landau levels and the quantum Hall effect which has been the subject of study in 
several references [e.g. ,2-4]. It is shown that w(a) satisfy the following relation: 

w(a)w(b) = exp(i—B.(a x b))w(a + 6), (4) 

where axii = 01^2 — 0261, which is exactly the q-commutation relation: 

w(a)w(b) = qw(b)w(a), (5) 
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where q = exp(i|S.(a x &)). We can write g in the form q = e n where 
-f = 2n-^. $ and <j> are the magnetic flux through the unit cell and the quan- 
tum of the flux respectively : 



$ = B.(axa$ = -. (6) 

e 

Here a and b are two perpendicular unit vectors which build the lattice unit 
cell, | a | = | b |= 1. Using (4) and (5) it is easy to show that the following 
combinations of the magnetic translation operators : 

= w w + w ® 

q-q- 1 

w(-a) + w(-b) 

3- — ! (°J 

q-q 

and g J ' 3 = w(b — a) = j z form the su q {2) algebra. The su q {2) algebra is the 
q-deformation of the lie algebra su{2) [5-8], with generators satisfying the com- 
mutation relations: 

[j 3 ,j±] = ±j±. (9) 

= Vf^- (10) 

We observe that (4), (5) and the su q (2) algebra are invariant under following 
transformations: 

w (a) -> 9(l)w(a), (11) 

w(-a) ^ g- 1 (q)w(-a), (12) 

or equivalcntly: 

j+^9(q)j+, (13) 

3-^g- 1 {q)3-, (14) 

3z->3z, (15) 
where g(q) is an arbitrary function of q. 

2. Areas distribution of closed random walks. 

Bcllissard ct al [9] derived the distribution of areas of closed random walks 
using noncommutative geometry. Their method was based on the Harper model. 
The Harper model Hamiltonian [10] is given by: 

H = ^2 w ( s ) = w ( s ) + w (b) + w(-a) + w(-b). (16) 
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We define the trace per unit area such that 

T(w(m 1 )w(rn 2 )) = 5mi+m 2 ,o- (17) 
Using (16) and (17) it is easy to see that : 

T(H N ) = Y, e ^ A(r) > ( 18 ) 
r 

where the sum is over the set of closed paths starting at the origin of length N, 
and A(T) is the algebraic area closed by T. If Sljv is the number of such closed 
paths, we have [9]: 

E" ^(|)exp(^) = ^ 1 E ex p(^) = ^ lT ^ JV )iT=t' ("J 

A— A max r 

where A max is the maximum algebraic area closed by T with nonvanishing prob- 
ability. Our aim is to find the formula for Pn(cl) where a is renormalized area 
a = ^, and to investigate its behaviour. 
From (7), (8) and (16) we obtain: 

H = w(a) + w(b) + w(-a) + w(-b) = 2i(q - q~ 1 )j y . (20) 

Thus we observe that the Harper Hamiltonian operator is a generator of the 
su q {2) algebra. If the eigenvalue of 73 is m, the eigenvalue of j z will be ±q m . 
Using (17) or its equivalent : 

T(w(c)) = 5 cfi . (21) 

one can show that: 

T(w(b- a)) = T£>(a)) = T(j y ) = 0, (22) 

a 

This result also follows from (9) and (10) directly. A straightforward calculation 
show that the eigenvalues of H have exponential form(see for example equ.(18)). 
To obtain the eigenvalues of if, we note that: 

1. E depends on the quantum number m, and is a function of q-deformation 
parameter 7(actualy 7 and N are not independent, 7 — j^, see equ.(19), and E 
can be a function of both of them) . It is also has exponential form ( for more 
details see appendix A) and satisfies in equ.(22). 
From mathematical point of view, if a series of functions 

00 

E-W- t )- ( 23 ) 

n=l 
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is convergent for any x in some interval, then this series defines a function of 
variable x, f(x) = J2^=i fn(%) m this interval. The importance of a series (e.g. 
a solution of some differential equations) lies in its convergence. For example, 
if I be an integer then the series expansions of the solution of the Legendre 
equation Pe(x) and Qe(x) are convergent, and they convert to polynomials. 
The same statement is also true for the hypergeometric functions and the other 
polynomials. 

We know that the energy of a physical system must be finite, therefore the 
series (6) in [9] is convergent and is a function of 7 (or N) which we will try to 
determine it. 

On the other hand we know that the operators j y and j z of su(2) algebra do not 
commute and can not have common eigenvectors, but their eigenvalues are the 
same. After deformation the eigenvalues of the j y and j z will be not the same, 
but still their eigenvalues will be discrete labelled by integer quantum number 
m. 

2. We know the eigenvalues of H, in the limit of large N [9]. 
Therefore we can write : 

E m = ie""»W. (24) 
If we choose u m (j) = —rrij + h(-f), then we have: 

E m = ±f{ 1 )e- m l, (25) 
where /(j) = cx.p(h(j)). As the candidate for 7(7) we take : 

f( 1 ) = 2{e { ™&& ) ' +e~ i{ ^^ ) ). (26) 

We show that its N — > oo is fully consistent with the result of [9] . We tested it 
numerically (see below). Then we will have: 

g ro = ±4cos( . JLj e- 2 * 1 , (27) 
smh(^) 

By introducing the new variable I 

m = 2t+l, (28) 

We have: 

TP. \-A nr.c( 

■ Nf 

■ 4 

Expanding Ei for small values of 7, 7 — > 0, we obtain: 



(2£+l)7 



Ei = ±4 cos( . -*- ) e --*-. (29) 
sinh(-ji) 



E t = ±4[(1 -J^)-(2t+l)l + [{21 + lf + l]t- 0(7^)], (30) 
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The term -^h? can be omitted, because it tends to zero when TV — > oo ( 7 — > 0), 
and equ.(30) reproduces the eigenvalues used by Bellissard[9]. Now T{H N ) is 
given by: 

tx^hEE^))^)' ( 31 ) 

where 7 = j| is the multiplicity per unit area. This leads to: 

aN+1 x x 

T ^^idwr eK( ^kr) > - (32) 

for even N. This gives us the exact characteristic function for the distribution 
of area. In the limit of large N we obtain: 

4N+I x , (x\2 1 

which is also derived by Bellissard [9]. From (19) and (32) we have: 

Q N = T(H N ) | 7=0 = — cos w (-), (34) 

This gives the total number of closed paths of length N. In the limit of large N 
we have: 

4W+1 1 
VN = ^ + (Ojf)). (35) 

Using (19), (32), (34) and the normalization condition we can obtain a series 
expansion for the probability distribution: 

^N\a) — ^rjv Z^ m =o 2^k=o 2^e=o 



{ 1)&m2 )[2e ) ((2m + 2fc + l)2 + 16 a 2)2 m+ 2 i ^ TO + 2 )' W 



We have plotted Pjv(d) as given by equ.(36) against exact enumeration [11] in 
figl. In the summation of equ.(36) we have allowed m to range from to 2 
and k from to 7000 for fig.l.(a) and from to 60000 for fig.l.(b) and (c). 
Comparing fig.l.(a) and (b), we find that a better result will be obtain if we 
increase the range of summation over k. The same statement is also true for m. 
m = and m — 1 give the equ.(9) in [9] and the correction term made by 
them respectively. Some of the coefficients are presented in Appendix B. 
In conclusion we observe that the expression presented for the spectrum of H 
using quantum groups is in good agreement with the results for N — > oo and 
also with exact enumeration. If there is other candidate for f(j), which full fills 
all the requirements, then the difference between it and our expression will be 
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negligibly small. 
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Appendix A. The Harper Hamiltonian is: 

H = e ai + e a2 + e- ai + e~ a2 = J2 ^ . ( 37 ) 

where [a^a,-] = 1^(7), e a *e a ^ = e ai+a ->' e^K^l and T(e a ») = S atfi . then: 
therefore we have: 

T(H N ) = ^e Um(7) . (39) 

m 

u m (-f) must be determined using boundary conditions. 

Appendix B. The first four coefficients are as follows: 

G N = 1 (40) 

Gf = - — . (41) 
C N - 1 1 + 1N ~ 1 (A 2 ) 

c n_ 11 1^-1 1 (N-l)(N-2) 
3 720 7V 5 48 N 5 48 N 5 1 ' 



Figure captions. 

Fig.l. The exact and approximate distribution for areas of closed random paths. 
The exact distribution is plotted for iV = 16. The solid points were obtained 
through exact enumerations [11]. 
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